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Abstract—This paper presents a review of complex 
coordinate approaches used in electromagnetics, which 
are based on extending the real coordinate space to a 
complex coordinate space. After a general discussion of 
methods used in the literature, two special problems are 
discussed in detail: non-Maxwellian perfectly matched 
layer (PML) designs in the context of finite element 
method, and perfectly matched double negative (DNG) 
layer design. These are supported by some numerical 
simulations. 


Index Terms—Complex coordinate, complex space, 
complex transformation, complex stretching, complex 
point-source, complex beam, complex ray, perfectly 
matched layer, double negative layer, metamaterial, 
transformation electromagnetics, transformation optics, 
anisotropic medium, transformation medium, 
electromagnetic scattering and radiation, computational 
electromagnetics, finite element method. 


I. INTRODUCTION 


Students or new beginners in engineering or physics 
usually come up with questions like: Are complex numbers 
related to physical reality? How are they used in real life 
applications? These are indeed deep questions with probably 
no trivial answers. It is usually believed that real numbers 
have physical reality because they are attributed to physical 
measurable quantities, but complex numbers have auxiliary 
or artificial nature in physics and can best be described by 
formal mathematics. However, a number (either real or 
complex) has a physical meaning if it can be linked to a 
theory in physics or engineering. Complex numbers, just like 
real numbers, are elements of mathematics, and can be used 
to model physical phenomena. They have a wide range of 
applications, such as signal representations, Fourier analysis, 
phasor analysis, etc. in engineering and physics. 

The purpose of this review paper is to present a unifying 
point of view that discusses the use of complex coordinates 
or complex space in electromagnetics. The key role of a 
complex space point is to introduce loss (dissipation) or gain 
into the problem of interest. This property can be used 
directly, for example to model losses in a medium, or can be 
utilized to create new functionalities with different purposes. 


The first use of the concepts of complex point-source (i.e., 
point source located at a complex space point) and complex 
rays can be traced back to Keller’s 1958 paper [1], followed 
by similar approaches to solve various scattering and 
diffraction problems in the literature [2]-[17]. It was shown 
that the concept of complex point-source can be used to 
generate Gaussian beams [3], [4]. The Green’s function 
expression e~/*" /r evaluated at a point P, where r is the 
distance between the source located at a complex point and 
the observation point P, represents the field of a Gaussian 
beam. Beams that are created by complex point-sources are 
exact directional fields determined by an analytical extension 
of the Green’s function into the complex space. In other 
words, the main mission of shifting points to complex space 
is to enable point sources to radiate in specific angles by 
providing decay in other angles. The field patterns of 
complex-point sources are similar to real-point sources 
except for the exponential angular decay. Note that the decay 
feature associated with complex coordinates is utilized to 
model sources with directionality. But, this decay feature 
was also used to solve scattering problems in lossy media by 
using the concept of complex ray tracing [5], [6]. Apart from 
these studies, Boag and Mittra showed that the matrix size of 
the method of moments solution of electromagnetic 
scattering problems can be reduced by representing the 
scattered field in terms of a series of beams produced by 
multipole sources located in complex space [16], [17]. The 
directional nature of the fields radiated by multipole sources, 
located in complex space, ensures that the coupling between 
distant parts of an object is low, and hence, matrix elements 
corresponding to low coupling can be eliminated. 

Perhaps, the most widely used application of complex 
coordinate transformation is the design of perfectly matched 
layer (PML), which is used as a numerical absorber in mesh 
truncation of finite methods (such as finite element method 
(FEM), and finite difference time domain (FDTD) method) 
for solving open-region radiation and scattering problems. 
Such a complex coordinate transformation is known as 
complex coordinate stretching in PML nomenclature. With 
the help of complex coordinates in the PML region, outgoing 
waves are replaced by exponentially decaying waves so that 
outgoing waves are absorbed without any reflection 
irrespective of their frequency and angle of incidence. The 
main superiority of the PML approach over absorbing 
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boundary conditions is that PMLs have the ability to 
minimize white space by locating it to the close proximity to 
the surface of the volume of scatterers or sources. In fact, the 
PML concept was implemented by means of different 
approaches and interpretations in the literature. It was first 
introduced by Berenger [18] in the context of the FDTD 
method, by using a split-field formulation of Maxwell’s 
equations in Cartesian coordinates. This approach yields 
non-Maxwellian fields within PML. Afterwards, Sacks et al. 
realized the Maxwellian PML in Cartesian coordinates as an 
anisotropic layer with permittivity and permeability tensors 
in the context of FEM [19]. Gedney used the concept of 
anisotropic PML medium in FDTD [20]. These anisotropic 
PML approaches were originally introduced in Cartesian 
coordinates. Kuzuoglu and Mittra extended this concept to 
cylindrical and spherical coordinates [21], and to the design 
of conformal PMLs using a local curvilinear coordinate 
system [22], in the context of the FEM. It was shown by 
Chew et. al. that the above approaches are equivalent to a 
more general approach, called complex coordinate stretching 
[23], [24]. This non-Maxwellian approach is based on the 
analytical continuation of the field variables to complex 
space. Teixeira and Chew showed that anisotropic PML can 
be realized by applying a complex coordinate stretching 
through the mapping of the non-Maxwellian fields obtained 
during the complex stretching to a set of Maxwellian fields 
in an anisotropic medium representing the PML [25]-[27]. 

The above-mentioned PML implementations, especially in 
the context of FEM, are based on the design of an anisotropic 
medium having suitably defined constitutive parameters. In 
addition, all of them were proposed in rectangular, circular 
or spherical domains that do not have arbitrary curvature 
discontinuities. Although the concept of local coordinates 
was used for designing conformal PMLs, its implementation 
is difficult especially for arbitrary geometries. Ozgun and 
Kuzuoglu proposed an approach, called locally-conformal 
PML, based on complex coordinate stretching to design 
conformal PMLs over challenging geometries, especially 
having some intersection regions or abrupt changes in 
curvature, in an easier and flexible manner [28]-[30]. Such 
conformal PML domains are vital to decrease white space, 
thereby to decrease computational load. This method does 
not employ the concept of anisotropic medium, and only 
replaces the real coordinates with their complex counterparts 
obtained by a complex coordinate transformation. Since the 
transformation is locally defined, it has no explicit 
dependence on the differential geometric characteristics of 
the PML-free space interface, and hence it becomes possible 
to model arbitrary PML geometries. Another method, called 
multi-center PML method, was also proposed for the design 
of conformal PMLs [31], [30]. Both methods will be 
summarized in Sec. II. Due to its simple and flexible nature, 
the locally-conformal PML method was also employed to 
divide a computational domain into a number of subdomains 
in domain decomposition methods to eliminate artificial 
reflections arising from the abrupt division of the domain 
(32]-[35]. It was also used to truncate infinitely-long 
conducting structures in a finite-sized domain to be able to 
achieve FEM modeling of high frequency diffraction 
problems [36]. 


In recent years, the coordinate transformation approach 
was also employed in designing metamaterial structures. In 
[37], Kuzuoglu extended the concept of a double negative 
(DNG) medium (an isotropic medium with negative 
permittivity and permeability), which was proposed by 
Veselago [38], to a perfectly matched DNG medium by 
means of complex coordinate transformations. He 
demonstrated the anisotropic generalizations of the Veselago 
medium. This approach will be explained in Sec. III. The 
concept of coordinate transformations was also combined 
with the transformation optics (TO) (or transformation 
electromagnetics (TEM)) approach. The TO approach is an 
intuitive and systematic approach to the design of 
metamaterials (or application-oriented transformation 
media), and it is based on “real” coordinate transformations 
to manipulate the behavior of electromagnetic waves in a 
desired manner. It uses the form invariance property of 
Maxwell’s equations under coordinate transformations [39]. 
When the spatial domain of a medium is modified by 
employing a coordinate transformation, this medium 
equivalently turns into an anisotropic medium in which 
Maxwell’s equations have the same mathematical form. The 
constitutive parameters of the anisotropic medium are 
obtained by the Jacobian of the coordinate transformation. 
The central concept of TO can be traced back to Bateman’s 
study in 1910 [40]. It was also discussed in the context of 
finite methods in [41], [42]. Over the last decade, the TO 
approach—which is based on real coordinates— was 
employed to the design of various optical and 
electromagnetic structures, such as invisibility cloak [43]- 
[45], electromagnetic reshaper [46]-[49], waveguide 
miniaturization [50]-[52], lenses [47], [53], [54], 
concentrators [55], [56], holes [57], antennas [58], [59], etc. 
Other than these applications that aim to design devices with 
new capabilities, the TO approach was also used to design 
computational materials (called software metamaterials) to 
overcome some numerical difficulties that arise in 
computational electromagnetics problems, such as modeling 
large-scale or multi-scale electromagnetic scattering 
problems [60], [61]; modeling curved geometries that do not 
conform to a Cartesian grid [62], [63], and modeling 
stochastic electromagnetic problems modeling uncertainties 
[64]-[67]. In recent years, the TO approach was extended, 
via complex extension of the spatial coordinates, to achieve 
field-amplitude control [68], to create single-negative [69] 
and PT-symmetric [70] transformation media, to handle 
non-Hermitian optical problems [71], and to design lenses 
[72]. Also, the use of complex coordinates was investigated 
to understand the effect of inhomogeneous media on wave 
propagation [73]. In these applications, the transformation 
medium involves both dissipation and gain due to material 
parameters that are complex-valued functions of position. 

The organization is as follows: In Sec. I, locally- 
conformal and multi-center PML approaches are 
summarized. In Sec. II, the design of perfectly matched 
DNG medium is presented. In Sec. IV, some conclusions are 
drawn. Throughout the paper, the suppressed time 
dependence of the form exp(jwt) is assumed. 
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II. PERFECTLY MATCHED LAYER (PML) DESIGN 


This section presents two PML approaches based on 
complex coordinate approaches: locally-conformal PML and 
multi-center PML. These are non-Maxwellian PML methods 
that do not require any artificial anisotropic medium and 
coordinate system. These methods use specially- and locally- 
defined complex coordinate stretching, and hence, make 
possible the easy design of conformal PMLs having 
challenging geometries. 

Before the details of these special methods, it is useful to 
discuss the principal idea in a classical PML approach based 
on complex coordinates. Consider first the simplest case as 
shown in Fig. 1, where a TM, plane wave (ie., 
E=4,E,(x,y)) with arbitrary frequency and direction of 


propagation is incident to an interface (T ) that separates two 
half spaces (x >0 and x<0). The PML and free space (FS) 


Qi ={(x,y) |x > of 


Q = {(x, y) |x < o} , respectively. In Qrs, the plane wave is 


regions are defined as and 


given by 
Elayne, (1) 
where k, =kcos@, k, =ksiné, k is the wave number and 


Ois the angle between the direction of propagation and the x- 
axis. The purpose of the PML approach is to transmit the 
plane wave into the PML region Qpm without any reflection, 
and to attenuate it in the +x direction. Therefore, it is 
stipulated that the plane wave in Qpm is given by 


E (x,y) = {Qe e ; 2) 
where f(x) is a function satisfying the following two 
conditions: 

e JOH 


e f(x) decreases monotonically for x >0. 


For example, f(x) can be chosen as 


f (x) = a rcot z (3) 


where œ >0 is a constant. Note that the angular dependence 
of the decay characteristic is needed to satisfy reciprocity in 
PML. Now, define a parameter a as 


a 
a=le-—, 4 
7k (4) 


By substituting (4) into (3), and then (3) into (2), we get 
E, (x,y) =e Tat IY (5) 

At this point, a complex coordinate x can be defined as 
xX=ax, (6) 


or equivalently, 


gage eg, (7) 





Fig. 1. 
interface. 


A TM, plane wave is incident on a PML-free space 


Note that, the real coordinate is transformed to a complex 
coordinate by just adding an imaginary part to the real 
coordinate x. With this complex coordinate, (5) becomes 


E, (,y)=e™e™". (8) 
This equation satisfies Helmholtz equation 
OE, OE 


g t p ET. (9) 





everywhere on R? , provided that 


Xx, x<0 


x= ; 10 
x a s50 (10) 


One attractive consequence of this interpretation is that it 
allows to solve the standard Helmholtz equation with a 
standard finite element method except that the real 
coordinates are simply replaced by complex counterparts. 
The above formulation is based on the fact that the PML 
region is a half space. Indeed, the PML region must be a 
bounded region enclosing all sources and objects. The above 
formulation can easily be extended to a rectangular layer in 
2-D Cartesian coordinates, a part of which is shown in Fig. 
2(a). The inner and outer boundaries are represented by Tj 
and I, respectively. A successful PML design ensures that 
the waves incident to Tj from the free space domain Qgs are 
transmitted into Qpmi without any reflection in order not to 
affect the fields within the free space region, and they decay 
as they traverse the PML region toward the outer boundary. 
Since the magnitude of the field reduces to a negligible value 
at the outer boundary, this boundary can be modeled as a 
perfect electric or magnetic conductor. In Fig. 2(a), the PML 
region is represented as a union of three regions. In regions I 
and II, the waves must be attenuated along the x- and y- 
directions, respectively. However, in the corner region III, 
the waves must decay in both directions. Hence, both 
coordinates are transformed to complex coordinates to 
attenuate the waves in an oblique direction that is a 
combination of the x- and y-directions. The complex 
coordinate transformation can be defined as follows: 
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Opm 


(a) (b) 


Fig. 2. Classical PML methods: (a) Cartesian, (b) Cylindrical. 


Region I: i=x+“(x-%) and j=y (11.a) 
J 
s 7 a, 
Region I: *=x and y=y+—(y-y,) (11.b) 


Jk 


Region II: ¥=x+(x-x,) and = y+2(y-y,) (11.0) 
Jk Jk 
Note that œ, and œ, might be chosen as identical if the rate 


of decay in both directions is desired to be the same. With 
this transformation, the plane wave expression is given by 


— jk,ax „—jkyy 
ere Š 


e EA e jkyby 


in region I 


in region II , (12) 


, in region III 


where a=1+(a,,/ jk) and b=1+(a@,/ jk). 

A similar approach can be developed if the PML region 
surrounds a circular spatial domain as shown in Fig. 2(b). 
Using the cylindrical coordinates, the waves are attenuated in 
the radial direction by defining the transformation as follows: 
(13) 


= a 
r=rt K (r-n). 
Note that these classical PML approaches fail if the PML 
domain is of arbitrary shape, because inner and outer PML 
surfaces must be defined on “constant coordinate surfaces”. 
The locally-conformal and multi-center PML methods 
(described below) overcome this difficulty, and the inner and 
outer PML surfaces need not be identical, providing a great 
flexibility in designing arbitrarily-shaped PML regions. 


A. Locally-conformal PML approach 


The method is illustrated in Fig. 3(a), which shows the 
spatial region occupied by a convex PML region (Qpmr). The 
PML can be designed as conformal to an arbitrary volume 
(Q) containing sources and objects. The PML is realized by 


mapping each point P in Q w, E R? to P in complex PML 
region Q pmu € C?, through the complex coordinate 


transformation as follows: 





Fig. 3. Conformal PML methods based on coordinate stretching: (a) 
Locally-conformal PML, (b) Multi-center PML. 


For+— f(E)é;, 


F (14) 


where re% and fFeC° are the position vectors of the 
points P in real space and P in complex space, respectively; 
k is the wavenumber; and € is the parameter defined by 


[pn 


where r, €T; is the position vector of the point P, located on 


, (5) 








the inner PML boundary Tj, which is obtained by solving the 
following minimization problem: 


min|[r—r, | . 


ner; 


(16) 


This implies that € is the shortest distance between the inner 
boundary and the corresponding PML point. Eq. (16) results 
in a unique r, because T; is the boundary of the convex set. 


Moreover, â z is the unit vector defined by 


a, =(r-r,)/¢, 


and f (é) is a monotonically increasing function of é as 


(17) 


follows: 


Forum for Electromagnetic Research Methods and Application Technologies (FERMAT) 


fe 


mir, -r, 


(18) 


m1 ? 








where gis a positive parameter (practically, 5k < œ <15k if 
PML thickness is between A/4 and A/2 where A is the 
wavelength), m is a positive integer (typically 2 or 3) related 
to the decay rate of the magnitude of the field inside PML, 
and r, is the position vector of the point Pp which is the 


intersection of the line passing through r and r,. Note that 
|r, —r,|| is the local PML thickness for the corresponding 
PML point. 


B. Multi-center PML approach 

This method (shown in Fig. 3(b)) is basically a 
generalization of the locally-conformal PML method. It is 
implemented by choosing a finite number of center points 
(such as P; and P2) inside the volume Q. The unit vectors (â, 


and â, ) emanates from these centers in the direction of the 
PML point, P. The unit vector a, , shown by the dotted line, 
is equivalent to a; in (17). Furthermore, 8; and 02 are the 


angles between â, and the unit vectors da, and 4d,, 


respectively. The complex coordinate transformation, which 


maps each point P inside the PML region to P in complex 
PML region, is defined as follows: 


m 
N = 
a X i [r Tai 
. i 
jk mir, -r 


a,i 








a, , (19) 


m-l `i 


r=r+ 








where Q@and m are the same as before, N is the number of the 
centers, and w; is the real weight value assigned to each 
center P;. The sum of all weights should add up to 1, and 
each weight w; is inversely proportional to the angle G, 
which is less than or equal to 45°. The value 45° is the 
threshold value in assigning the weights, and has been 
determined empirically by means of systematic numerical 
experiments. The weight selection scheme assigns the 
highest priority to the center whose unit vector â, makes the 


smallest angle with â,. To achieve smooth decay inside the 


PML region, the centers should be chosen in such a way that 
each PML point should have at least one center, whose angle 
@ is less than or equal to 45°. In some smooth geometries 
(such as a spherical or cubical shell), a single center-of-mass 
point can provide reliable results. The center selection 
scheme is a straightforward task depending on the geometry 
of the PML region. 

Both locally-conformal and multi-center PML 
transformations can be implemented by just simply replacing 
the real-valued node coordinates inside PML by their 
complex-valued counterparts calculated via transformations. 
Another way is to model the PML region as an anisotropic 
medium whose constitutive parameters are obtained by the 
form invariance property of Maxwell’s equations under 
coordinate transformations. This approach will be 
summarized in Sec. II.C for the sake of completeness. 


Let us now discuss how the transformations achieve field 
decay inside the PML region. Consider an outgoing wave in 
the neighborhood of an arbitrary point Pa, which can be 
locally represented as a superposition of plane waves. A 
typical representative is given as follows: 


E(r)=4,¢"", (20) 
where a, is the unit vector along the direction of incidence, 
and a, is the unit vector denoting the polarization. Let 
E(i) be the transformed version of E(r) to complex space 
given by the expression 

E(F)=4,e%* . (21) 


By using the coordinate transformations, the field inside 
the PML region becomes 








E(r)=E(r) g(r), (22) 
where 
g(r)= exp] -f (é) â, â; | , (local-conformal) (23a) 
N ar-f. m 
g(r)=exp -> w; eal a -@, |, (multi-center) (23b) 
i m| Ri T Ti 











This function satisfies the following two conditions: 


(i) g w), „ =1 . It assures the continuity of field at 


the interface to get rid of numerical reflections. 
(ii) It is a monotonically decreasing function away 
from the PML-free space interface. 

Consequently, the coordinate transformations satisfy the 
following three conditions that are essential for a successful 
PML design: 

(i) the outgoing wave in the neighborhood of T; is 
transmitted into QpmL without any reflection, 

(ii) the transmitted wave decays monotonically within 
Qpm along the direction of the unit vector and away 
from the PML-free space interface, 

(iii) the magnitude of the transmitted wave is negligible 
on To. 

The performances of the PML methods are tested by FEM 
simulations considering the TM, radiation of a line-source 
inside an infinitely-long cylindrical region of hexagonal 
shape. The position of the line-source is assumed to be a 
“random variable” uniformly distributed inside the region, 
based on the assumption that the performance of the PML 
may depend on the position of the source. For this purpose, 
we utilize the Monte Carlo simulation technique, which is a 
method based on the use of repeated random simulations to 
extract the statistical parameters of stochastic problems, to 
show that the PML methods provide reliable and robust 
results irrespective of the source position. We randomly 
determine 2000 different source positions, and run the FEM 
program 2000 times. For each run, we calculate the mean- 
square error by comparing the electric field calculated by the 
FEM, which is truncated by either locally-conformal or 
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multi-center PML, with the analytical solution. In Fig. 4, we 
plot the error scatter plot that shows the error value at each 
source position, as well as the error histogram. We also show 
some Statistical data (i.e., mean, variance, etc.) on the plots. 


C. Maxwellian PML implementation via form invariance 
property of Maxwell’s equations 


The non-Maxwellian PML methods described above are 
implemented by simply interchanging the real coordinates by 
complex-valued coordinates in a standard FEM code. They 
do not model the PML as an anisotropic medium. However, 
the approaches are indeed equivalent to anisotropic medium 
design through the use of form invariance property of 
Maxwell’s equations. The medium in which the 
transformation is applied turns into an anisotropic medium 
with constitutive parameters obtained by the Jacobian of the 
transformation. Transformed fields in the anisotropic 
medium satisfy the original form of Maxwell’s equations. 
This approach is briefly summarized below. 

If the original medium before the transformation is an 
isotropic medium with parameters (€, 4), then the 
permittivity and permeability tensors of the anisotropic PML 
medium are determined by 


Fay, =EA and py, =HA, (24) 


where 

A=6, Ty", (25) 
and J= O(%, J, Z)/0(x, y, z). Here, J is the Jacobian tensor 
and 6; is its determinant. If the original medium is an 
anisotropic medium with parameters (&, 7’), then the 
parameters of the anisotropic PML medium are obtained by 


Eng =6,(T") -B-(T*), 
Ta aT] AT 


Assuming that the wave propagation is along the infinite z- 
axis (i.e., no z-variation), z-dependent off-diagonal terms of 
the Jacobian tensor become zero (i.e., 
dx/dz = 09 /dz =0Z/dx =0Z/dy=0), as well as the z- 


dependent diagonal term becomes unity 0Z/dz =1. 


(26a) 


(26b) 


Based on the form invariance property of Maxwell's 
equations, the electric field satisfies the following 3-D vector 
wave equation in transformed and original coordinates, 
respectively, as follows: 


: (27a) 
: (27b) 


where V =[J"]"-V is the del operator in transformed space, 
E(f) is the original field in transformed coordinates, and 


E(r) is the transformed field in original coordinates. These 


fields are interrelated through the following principle of field 
equivalence: 
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Fig. 4. Monte Carlo simulations to test the performance of the PML 
methods: (Top) Mesh structure, (Middle) Error scatter plots, 
(Bottom) Histogram plots. 


E(r) = T -E(Č), 
H(r) =J" -H(@). 


(28a) 
(28b) 


For TM, case, the vector wave equations in (27a) and 
(27b) reduce to scalar Helmholtz equations in transformed 
and original coordinates, respectively, as follows: 


(29a) 
(29b) 


VE, +k’E, =0, 
V-(A,VE,)+r6A,,E, =0, 


where A, is the 2x2 tensor corresponding to the transverse 
part of A, ô, is its determinant, and A, is the (3,3) 


component of A. Because of the special form of the 
Jacobian tensor in TM, case, (28a) is expressed as 


E(%, 53) =E,(x, y). 


II. PERFECTLY MATCHED DOUBLE NEGATIVE LAYER 
(DNG) DESIGN 


In a double negative (DNG) medium, which is also called the 
Veselago medium, the wave vector k and Poynting’s vector 
P are anti-parallel, that causes negative refraction and 
reversal of Doppler shift and Cherenkov radiation. In this 
section, we model a perfectly matched DNG medium by 
using complex coordinate transformations, and show how to 
obtain anisotropic generalizations of the Veselago medium. 
It should be noted that although the coordinate 
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transformations employed in this section are eventually 
reduced to real-valued mappings, the ideas are very similar 
to those in the previous section. 


Let us define a complex coordinate mapping T :R* > C?, 
which applies the transformation only in x-direction, as: 


F=(x,y,z) > F=T(F)=(%5,D, (30) 


where 


x x<0 


X= x d y= d Z= 
Xx f s, dx an y=y and z =z 
0 


x>0 an 


where the stretching parameter s, is in the complex form 
s (x) = S(x) + js (x) . Hence, in the half space x>0, the 


nabla operator V in the complex space and the tensor A in 
(25) , respectively, become 





VWs, 0 0 
V=| 0 olv (32) 
0 0 1 
_ ji/s, O 0 
A=| 0 s, 0 (33) 
0 O s, 
because 
Dai T oat Ca (34) 
ox “o dy oy oz a 


Now, let (E,H) be the transformed fields (i.e., mapped 


versions of (E,H) to C*) by restricting them to rec, 
defined as the manifold 


P={(%, 5, C'|T(x, y,2)= (45,0, (ye R}. Thus, 


which is 


the governing equations in the region 
Q, ={(x, y, eR |x > o} become 

VxE=—jau,H (35.a) 

VxH = jae,E (35.b) 

Hence, the field distribution (E“,H*) inside the 

anisotropic medium occupying the region Q, satisfy 


Maxwell’s equations with constitutive parameters 7 = ,A 


and Z = E&A , and can be expressed as 


(36.a) 
(36.b) 
The stretching variable s, is in general a function of the 


space variable x. Let us assume that s, is a complex 
constant given by 


$,=a-jo (37) 


(Ue i 8 Bur ET 
Pater Te 


Frée*spaee 


k2 


y (ind) 


y (in d) 





5 0 5 
x (in 2) 


(b) 
Fig. 5. Transmission of a plane wave into the perfectly matched 
DNG medium: (a) Coordinate stretching in x-direction, (b) 


Coordinate stretching in x- and y-directions (Overlapping perfectly 
matched DNG layers) 


where œ and o 20 are real constants. In the region ©, , the 
space variable in the transformed space becomes 


¥ =(a-jo)x (38) 


In the limit as o +0, the anisotropic medium becomes 
isotropic. If @=1, the medium is free space (€=€) and 
H=Hp), whereas if @=-1, it is a Veselago medium ( 
€=—-€) and “=—fly). 

Let us now illustrate how this transformation creates 
negative refraction in the Veselago medium with €=—€) 


and £4=—{/,. We assume that a TM, plane wave (E,H) 
exists in the half-space region Q_ ={(x,y,z)e R |x <O}. 


Defining the propagation vector as 
k, =k (4, cos@ +4, sin 8) , the fields can be expressed as 


follows: 
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E=4, exp| — jk (xcos0+ ysinð) | (39a) 
H= Ifa, sinĝ-â, cos 6 |exp| — jk (xcos 6+ ysin 0) | 


i (39b) 


k, XE 
a (hE) 


where 7) is the intrinsic impedance of free space. The triplet 
(E,H, k,) is right-handed, as illustrated in Fig. 5(a). In the 


half-space Q,, the analytic continuation of E and H to 
complex space becomes 


fh 
ll 

& 
oO 


xp| — jk (Xcos 6+ fsin 8) | 
(40a) 
pl -jk ((—x)cos 6+ y sin a) | 


H= “fa, sin 8-4, cos Ø |exp| -jk (žcos6+ ysin 0) | 


Th 
= ala sin- â, cos 9 |exp| -jk ((~x)cos 0 + ysin 6) | 

Tho 
(40b) 


By using (36), the fields E° and H* are obtained as 


E“ =â, exp| — jk ((-x)cos 8 + ysin 0) | (41a) 

H= adia, sinĝ +â, cos 6) exp| —jk ((—x)cos 0+ ysin a) | 
h 

(41b) 


Defining the propagation vector as 


k, =k (-cos 6a, +sin 64, ) , these fields can be expressed as 


E* = 4, exp(— jk, -r) (42a) 


Ba rd 
H ae xk, ) 


(42b) 


We can conclude from the field expressions that the 
perfectly matched Veselago medium yields negative 
refraction and the left-handed triplet (E*,H*,k,), implying 
that the wave vector k, and Poynting’s vector P are anti- 
parallel. 

So far, we assumed that stretching is applied in x 
coordinate only. This approach can be generalized by 
defining the mapping in all three coordinate variables. We 
define the coordinate transformation as follows: 


x x<0 

a [iss (x)dx’ x>0 cee) 
y y<0 

Ve f s (y)dy y>0 re 
Z z<0 

(43c) 


Z= z P s 
[s-a z>0 
0 


where s, =s, =s,=—Q— jo. Defining a new set of fields 


y 
E‘ =(-a- jo)E and H’ =(-a-— jo)H, this set of fields 
satisfy Maxwell’s equations in an anisotropic medium with 


T = h^ and Z= EA , where 


SS, /8, 0 0 


y 


> Il 
I 


0 5,5, / 8, 0 (44) 
0 


0 5,8, Is, 


Letting @=1 and o 40, we obtain A=+I depending 
on the octant where the point (x, y, z) is situated. 

Consider the 2D case where s, =1, as shown in Fig. 5(b). 
Let the region Q, ={(x, ye R? |x <0, y< o} be free space 


and assume a field variation in the form of a plane wave 
given by 


E =4, exp(— jk, -r) (45a) 
H=- (k,xE) (45b) 
Olh 


where k, =k (4, cosĝ +â, sin 0) . After extending (E,H) to 
complex space and obtaining (E,H), the fields in the DNG 


media are obtained via (E“,H*) —identified as (E,H) — as 
follows: 


(i) In the region Q, ={(x, yje R’ |x>0,y< o}, s,=-l 
and Sy =1. Then, 
A=-I (46a) 
E =â, exp} jk, r) (46b) 
H =—_(Exk,) (460) 
alh 


where k, =k (-cos Oâ, +sin 6â,) ; 
(ii) In the region Q, ={(x, ye R’|x<0,y> o}, s, =1 and 
Sy =-1. Then, 


R= (47a) 

E =â, exp(- jk, r) (47b) 

H=——(Exk,) (470) 
Olly 


where k, =k (â, cosĝ-åâ, sin 2) : 
(iii) In the region Q, ={(x, ye R |x>0,y> o} , s, =-1 


and sy = —1. Then, 


A=1 (48a) 
E=4, exp(—jk,-r) (48b) 
H=— (k, xE) (48c) 


Dy 
where k, = k(-4, cos 8-4, sin o) ; 


The expressions (45-48) can be interpreted as follows: 
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(i) The constitutive parameters in the regions Q, and Q, are 
E and 4o. The region Q, is originally free space, and the 
region Q, corresponds to the overlapping of the two 


perfectly matched DNG layers. It is remarkable that (44) 
implies 

(a) If two perfectly matched DNG layers overlap, the 
medium is equivalent to free space (2D corner region). 

(b) If three perfectly matched DNG layers overlap, the 
medium is again DNG (3D corner region). 
(ii) Inregion Q,, k, =—k,. 
(iii) If we choose a closed surface S enclosing the corner, the 
net power flow through the surface is zero, as can be 
concluded from the Poynting’s vectors in the four quadrants. 

As a special application of the approach introduced in this 
section, we can construct a perfectly matched anisotropic 
Veselago slab. Let us consider a Veselago medium 


occupying the region Q, ={(x, y,. ze R? |x, <x< x}, and 
the region outside the slab is free space. The complex 
coordinates (x, ¥,Z) can be expressed as 


x x< xX, 
X=4 x, +5, (x—x,) X,<X<X,, 

X, +8, (X, —X,)+(X-X,) X>X, 
Y=y,Z=z (49) 


The stretching variable s, is chosen as s,=—-Q@-jo, 
where @>0. In the limit as o —> 0 , (49) becomes 


x x< xX 
pai gE] x <x<x, (50) 
X, — A(X, —X,)+(X—- ax) X>X, 


Note that when @ =1, the slab becomes perfectly matched 
isotropic Veselago slab with £ = —£) and 4 = -4 . 

The slab is perfectly matched to free space, and acts as a 
perfect lens, with more general focusing properties, under 
certain conditions. To illustrate the lens-like action, let us 
consider the 2D TM, case, where a line-source, extending to 
infinity along the z-direction, is represented by the current 
density as follows: 


J =4,16(r-r,) (51) 


where J is a constant (amp/m). The source is located at 
r,=(x,,y,), where x, <x. Starting with Maxwell’s 
equations, the scalar Helmholtz equation for the z-component 
of the radiated electric field ( E, ) can be written as 


VE, +kE, = jau,l6(r-r,) (52) 
The analytical solution can be computed by 
k 
E, (xy)=- e (ke-r) (53) 


where H is the Hankel function of the second kind of 


zero-th order, 77 is the intrinsic impedance of the medium, 


and |r—r,|= (x—x,)'+(y—y,) . The mapped version of 





E,(x,y) to complex space is E, (x,y), and using (50) the 


solution in the presence of the perfectly matched Veselago 
slab becomes 


X<X, 


nofe ao) 
kn 


E, (x,y) ot Ho [« (x, —A(x— x,)—x, f +(y- ny) 





nP (r fazat o | X> X 
(54) 


It is clear from (54) that a focus within the slab (i.e., 
Xi < Xp < x ) will occur at 


1 1 
X,, =(+—)x,-—x, (55) 
a a 
provided that the following inequalities are satisfied 
xX, <x, (56) 
xx, < Alx — x,) (57) 


The inequality (56) always holds since the source is 
outside the slab, and (57) indicates that the distance between 
the source and the slab must be less than the slab width ( 
d =x,- x) multiplied by @. Under these conditions, the 


position of the focus outside the slab will be 


X, =X,+(1+a)d (58) 


fo 


We can conclude that the distance from the dipole to the 
focus outside can be adjusted by changing the parameter @ 
in the anisotropic Veselago slab. In other words, the 
anisotropic slab rather than the isotropic Veselago slab 
provides flexibility in the lens design. To demonstrate the 
concepts in the design of anisotropic Veselago slab, we plot 
the contours representing the real part of the electric field in 
Fig. 6 for different levels of anisotropy, namely for various 
values of @. The parameters of the slab can be easily 
visualized on the graphs. 


IV. CONCLUSIONS 


We reviewed complex coordinate approaches used in 
electromagnetics. A complex coordinate transformation is 
based on extending the real coordinate space to a complex 
coordinate space. We focused on two complex space PML 
methods (locally-conformal and multi-center PML methods) 
which are used for mesh truncation. We concluded that these 
methods are implemented by just replacing the real 
coordinates by complex coordinates inside the PML, and 
they provide great flexibility to design conformal PMLs over 
arbitrary geometries. Afterwards, we discussed the design of 
perfectly matched double negative (DNG) layer using 
complex coordinates, and illustrated the focusing feature of 
the resulting metamaterials. The proposed designs were 
validated through some numerical simulations performed by 
the finite element method. 


Xp <X<X, 
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Fig. 6. Perfectly matched Veselago slab with d = 2A: (a) Line- 
source at xs = 11 in free-space (& = -1) ; (b) Isotropic DNG medium 
with @ = 1. Focus is at xf = 3A inside, xj = 5A outside; (c) 
Anisotropic DNG medium with a = 2. Focus is at x = 2.5A inside, 
Xfo = 7X outside; (d) Anisotropic DNG medium with @ = 3. Focus is 
at xfs = 2.33A inside, xfo = 9A outside. (A = 1m) 
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